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Abstract 

We consider a two-queue polling model with switch-over times and fc-limited service (serve at most ki 
customers during one visit period to queue i) in each queue. The major benefit of the fc-limited service 
discipline is that it - besides bounding the cycle time - effectuates prioritization by assigning different service 
limits to different queues. System performance is studied in the heavy-traffic regime, in which one of the 
queues becomes critically loaded with the other queue remaining stable. By using a singular-perturbation 
technique, we rigorously prove heavy-traffic limits for the joint queue-length distribution. Moreover, it is 
observed that an interchange exists among the first two moments in service and switch-over times such that 
the HT limits remain unchanged. Not only do the rigorously proven results readily carry over to N(> 2) 
queue polling systems, but one can also easily relax the distributional assumptions. The results and insights 
of this note prove their worth in the performance analysis of Wireless Personal Area Networks (WPAN) and 
mobile networks, where different users compete for access to the shared scarce resources. 
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1. Introduction 

This note studies a two-queue fc-limited polling 
model in heavy-traffic (HT), in which one of the 
queues becomes critically loaded with the other 
queue remaining stable. Under the fc-limited strat¬ 
egy the server continues working until either a pre¬ 
defined number of ki customers is served at queue i 
or until the queue becomes empty, whichever oc¬ 
curs first. The fc-limited policy is easy to imple¬ 
ment, augmenting the exhaustive strategy (where 
ki = oo) by a single “knob” for each queue that can 
tune performance to be efficient or fair. 

Significant gains in system performance have 
been observed in many application areas, e.g. com¬ 
munication S i) and production systems (see, 
e.g. HI), by implementing the fc-limited strat¬ 
egy instead of the traditional exhaustive and gated 
strategies (see, also, Q for a survey of applications 
of polling systems). Our specific interest for this 
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model is fueled by performance analysis of (Blue¬ 
tooth) Wireless Personal Area Networks (WPAN) 
and mobile networks, where different users compete 
for access to the shared scarce resources. In case of 
a Round-Robin type of scheduling, polling models 
with limited service policies naturally emerge (see 
i)- Proper operation of these systems and the 
scheduling therein is particularly critical when the 
systems are critically loaded. 

A mathematical challenge is that the fc-limited 
strategy does not satisfy the branching property for 
polling systems, which significantly increases the 
analytical complexity 11|. A testament for this 
statement is that - notwithstanding the wealth of 
literature on polling systems - the problem of deter¬ 
mining the steady-state distribution for our model 
is one that remains open. The statement is rein¬ 
forced as well by the fact that heavy-traffic lim¬ 
its have only been rigorously proven for branching- 
type service disciplines with very few exceptions. 
Non-branching policies differ from their branching 
counterparts both in the techniques useable and 
the resulting HT behavior. That is, the tech¬ 
niques used for branching-type policies typically 
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rely heavily on explicit steady-state results for sta¬ 
ble systems; however, corresponding results for non¬ 
branching policies are not available. Moreover, non¬ 
branching systems can possess both stable and in¬ 
stable queues which is very different from branching 
policies where all queues become instable simulta¬ 
neously leading to incomparable HT behavior. 

We rigorously prove HT asymptotics for two- 
queue fc-limited polling models with switch-over 
times via the singular-perturbation technique. 
That is, by increasing the arrival intensity at one 
of the queues until it becomes critically loaded, we 
derive HT limits for the joint queue-length distri¬ 
bution. The singular-perturbation technique can be 
typified by its intrinsic simplicity and its intuitively 
appealing derivation, although it requires some dis¬ 
tributional assumptions (see Q for a survey of ap¬ 
plications of the perturbation technique to queueing 
models). This technique enables us to present the 
first study rigorously proving HT asymptotics for 
polling models with a non-branching service strat¬ 
egy. It is proven that the scaled queue-length dis¬ 
tribution of the critically loaded queue is exponen¬ 
tially distributed, whereas the number of customers 
in the stable queue has the same distribution as the 
number of customers in a vacation system. More¬ 
over, both queue-length processes are independent 
in HT. The analysis can be directly extended to an 
Wqueue system {N > 2) with one queue becom¬ 
ing critically loaded, as well as one can relax the 
distributional assumptions. 

As a by-product, we observe a close similarity 
between fc-limited polling systems with and with¬ 
out switch-over times. This novel result, which 
is reminiscent of the variance absorption result in 
branching-type polling systems Q, shows that in 
HT the first two moments of the switch-over times 
can be effectively absorbed into the corresponding 
moments of the service times without impacting 
the HT asymptotics of the instable queue. In this 
regard, new fundamental insight is gained which 
should prove useful in the analysis of non-branching 
polling systems. 

The approach in the present note has its origin 
in [^, where systems with zero switch-over times 
are studied. At face value the extension to nonzero 
switch-over times may seem a small one, however 
this extension impels us to, considerably, modify 
and extend the analysis in [^. This reveals itself 
clearly in the determination of the parameter of 
the exponential distribution for the scaled queue 
length of the critically loaded queue. In case of 


zero switch-over times this parameter follows di¬ 
rectly from the HT limit of the total workload in a 
standard MjGjl model and the averaging princi¬ 
ple. When analyzing a system with nonzero switch¬ 
over times, service capacity is obviously lost due 
to the switching and, as a consequence, one can¬ 
not directly state that the total workload equals 
the MjGjl workload and the aforementioned in¬ 
terchange of switch-over and service time moments 
is required. Lastly, we note that in many practi¬ 
cal applications the switch-over times may be sub¬ 
stantial and that the presence of these switch-over 
times may be crucial for the operation of the system 
(see, e.g., [l^)- In the interest of space, we do not 
give all proofs in the current note, but we focus on 
the main theorems, the theoretical insights result¬ 
ing from these theorems and the observed similarity 
between polling systems with and without switch¬ 
over times. 

2. Model description 

We consider a polling model consisting of two 
queues, Qi and Q2, that are alternately visited by a 
single server. When ki customers have been served 
or Qi becomes empty, whichever occurs first, the 
server switches to the other queue. Then, an expo¬ 
nentially distributed switch-over time Si with pa¬ 
rameter <Ji is incurred. If the other queue turns out 
to be empty, the server switches back and serves, 
again, at most ki customers. The situation where 
both queues are empty will turn out to be negligible 
in the HT limit, which is considered in this note. 
A possible scenario is that the server waits until 
the first arrival and switches to the corresponding 
queue (say, Qj) to start another visit period of at 
most kj customers. Customers arrive at Qi accord¬ 
ing to a Poisson process with intensity A^. We as¬ 
sume that the service times of customers in Qi are 
independent and exponentially distributed with pa¬ 
rameter fii. We denote the load of the system by 
p = Pi + P 2 i where pi = \ijp,i. The utilization Ui 
for queue i, i = 1,2, is defined as follows 


where S is the total switch-over time in a cycle, i.e. 
E[S'] = l/cTi -|- 1 /ct 2. The cycle time G is defined 
as the time between two successive visit beginnings 
to Qi, with mean E[S']/(1 — p). A known result for 
this model (cf. 0) is that a necessary and sufficient 
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stability condition for queue i reads Ui < 1, i = 1,2 
Q. This condition can be rewritten as follows, 

A,E[C] < h. (2) 

In words, this means that for a stable system the 
average number of type-i customers arriving in a 
cycle is smaller than the service limit ki, i.e., the 
maximum number of type-* customers served in a 
cycle. Note that this system might have one sta¬ 
ble queue (with utilization less than 1) and one in¬ 
stable queue (with utilization greater than 1), in 
contrast to the commonly studied branching-type 
polling systems where all queues become instable 
simultaneously (see, for example. Resing [HI). 


3. Analysis 

We study the limiting behavior of the model in 
HT. In particular, we take the limit such that the 
utilization of Q 2 reaches its critical limit I, while 
Qi remains stable. In this note we take the limit 
by increasing the arrival rate of Q 2 until it reaches 
its critical value while keeping the other parameters 
fixed. Note that increasing X 2 will also increase 
Ml, the utilization at the other queue. In order to 
prevent Qi from becoming instable, we impose the 
following condition: 

Ai/fci < A2/fc2- 


to Qi- We refrain from giving all transition rates 
as they are trivial to determine. 

The next step, after having defined the states and 
transition rates, is to give all balance equations so 
we can apply a singular-pertubation to them. We 
select one balance equation to illustrate the tech¬ 
nique, 

(Ai-|-A2-l-cri)p(0,n2,fci-l-I) = A2p(0,n2-I,fci-|-I) 

ki 

+ h) + cr2^>(0,n2, fci + /C2 + 2), (3) 

h^l 


for 712 > I. This equation relates the in-flow and 
out-flow of state (0, 712 , fci + I), which is the state 
where Qi is empty, Q 2 has 712 customers, and the 
server is switching from Qi to Q 2 - One can leave 
this state through an arrival at any of the queues, or 
because the switch-over time has ended. This state 
can be entered from state (0 ,712 — I, -|-I), through 
an arrival at Q 2 ^ or through a departure from Qi 
when in any of the states (1, 712 , h) for h = I,..., fci, 
or from state (0,7i2, fci -I- ^2 + 2). In the last case 
the switch-over time S 2 has finished, but since Qi is 
empty, the server starts switching back to Q 2 again. 

We now introduce the following definitions 


M2 k2 


Var[R'] = 4 + ^^, (4) 
M2 ^2 


The results still hold when other (combinations of) 
parameters are modified such that M 2 11- See Sec¬ 
tion 0] for a short discussion on this topic. The case 
where both queues become instable simultaneously 
requires a different analysis leading to different re¬ 
sults and is beyond the scope of this note (see Sec¬ 
tion IH). 

3.1. Perturbed balance equations 

The model is a continuous-time Markov chain, 
and we start by describing its states {ni,n 2 ,h), 
where Tii and 712 denote the queue lengths. The pos¬ 
sible server states are numbered h = I,2,...,/ci -I- 
^2 -I- 2 and should be interpreted as follows. States 
I < M < fci indicate that the server is serving the 
h-th customer during its current visit to Qi. State 
h = indicates that the server is switching from 
Qi to Q 2 . Similarly, values ki-\-2 < h < ki-\-l-\-k 2 
indicate that the server is serving the h — ki — I-th 
customer at Q 2 . The final state, h = fci -|- /c 2 + 2, 
indicates that the server is switching from Q 2 back 


where Var[S'] = ^ . The random variable B '2 

will be given an interpretation in Section 3] 

Perturbation. We increase the arrival rate 
of Q 2 as follows to its critical value, 

A 2 = - du, w>0,0<5<I. (5) 

in [±* 2 ] 

Substituting ([HI) in o, taking 7 = 2, and letting 
(5 4-0 implies that M 2 t I and, thus, that Q 2 becomes 
instable. An appropriate value for the constant w 
is chosen at the end of this section. 

Let ^ = (5772, and 

M(ni,^/d,h) = (5(/)„,,/,(^,(5), (6) 

for 0 < f = 0{1), h = 1,2,..., ki-\-k2-\-2. Note that 
once uj is chosen, 6 and the scaled variable ^ will be 
uniquely defined. The next step is to substitute 
and (IH) in the balance equations. For Equation m 
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we obtain. 


We obtain the density Po(0 of the scaled 
number of customers in Q 2 by solving the 
differential equation (|S]), where we use that 

= 1 and that 


^=1 

Taking the Taylor series with respect to S yields a 
set of equations. 

3.2. Main results 

In this subsection, we derive the limiting joint 
distribution of Ni and SN 2 , the queue length at Qi 
and the scaled queue length at ( 52 , as 5 j, 0. The 
three results are obtained by equating the 0 ( 1 ), 
0 { 6 ), and 0{5^) terms in the set of equations from 
the previous subsection. These steps follow the line 
of reasoning introduced in Q and, therefore, the 
details are omitted in interest of space. 

We define L^^\z) as the steady-state queue- 
length PGF of a single-server multiple-vacation 
queue with parameters Ai and /ii, with fci-limited 
service. The vacation distribution in this system 
equals the convolutions of ^2 service-time distribu¬ 
tions in Q 2 and of the switch-over times Si and 
82 . We also define Po{0 as the scaled queue-length 
distribution of Q 2 , which is yet to be determined. 
Equating the 0(1) terms in the perturbed balance 
equations gives the following result. 

Theorem 1. 

00 ki-\-k2-\-2 

E E <,.( 0 ^”^ = ^'‘’H^)^o(e). ( 7 ) 

ni—0 ^—1 

For this vacation system, the PGF of the queue- 
length distribution is known (cf. 0 ). 

Subsequently, we determine the unknown PoiO 
in © by solving the equations for the first-order 
and second-order terms in the perturbed balance 
equations. This results in 

OoP^'(0 + CiP''iO + 02P'(6 = 0, (8) 


where the constants Cq 

, Cl, C 2 are defined as fol 

lows. 






+ 

1 

E[S'l 

(- 

- [E[P'2] _ 2E[P']2] 

Oo — 




2E[P']2 

Cl = 

1 

/ 

>>1 ^ 

1 , C2=UJ. 

E[P'] 

1- 



V 

k-lj 



/ypo(e)de = i. 

Theorem 2. 

-Po(C) = (9) 

with 

V = - / s -. (10) 

+ife (1 - ft) nB?] 

A natural choice for w is w = 1/E[i32]j lead¬ 
ing to the commonly used scaling lim^j^o ^^2 = 
lim„ 2 -f-i(l — U 2 )N 2 . We refer to Section 0^ for a 
discussion on this choice and alternatives. By ap¬ 
plying the multiclass distributional law of Bertsi- 
mas and Mourtzinou [i| it directly follows that the 
scaled waiting time at Q 2 also follows an exponen¬ 
tial distribution with parameter \ 2 rj. 

Finally, the above has the following immediate 
consequence for the joint (scaled) queue-length dis¬ 
tribution in HT. 

Theorem 3. For Xi/ki < \ 2 lk 2 and A 2 = (1 — 
Pi — 5 )/E[B' 2 \, we have: 

limP[A^i < ni,5N2 < C] = ^(«i) (l - , 

54'0 

(11) 

where 77 is given by cni), and C(-) is the cumula¬ 
tive probability distribution of the steady-state queue 
length of a single-server multiple-vacation queue 
with parameters Ai and pi, with ki-limited service. 
The vacation distribution in this system equals the 
convolutions of k 2 service-time distributions in Q 2 
and of the switch-over times Si and 82 - 

4. Discussion 

4-.1. Interpretation 

The main result m has the following intuitively 
appealing interpretation, which we derive heuristi- 
cally below: 

1. The number of customers in the stable queue 
has the same distribution as the number of cus¬ 
tomers in a fc-limited vacation system. 

2. The scaled number of customers in the criti¬ 
cally loaded queue is exponentially distributed 
with parameter rj. 
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We now let Ai —>• A* and A 2 —t A 2 for <5 4 , 0, with 


3. The number of customers in the stable queue 
and the (scaled) number of customers in the 
critically loaded queue are independent. 

Property [T] follows from the fact that if Q 2 is in 
HT, then exactly k 2 customers are served at this 
queue during each cycle. 

For Property m we take a look at an alternative 
A:-limited polling system with the same arrival rates, 
no switch-over times, the original service times at 
Qi and adjusted service times at Q2 with mean 
£[^ 2 ] and variance VarjS^]. Since in the original 
system the two switch-over periods and S 2 are 
always accompanied with precisely ^2 service times 
at Q 2 in HT, the first two moments of the amount 
of time the server is utilized - either due to serv¬ 
ing or switching - in the original and alternative 
system is obviously identical. In Q it is proven 
that the distribution of the scaled workload within 
the alternative system without switch-over times - 
and, thus, the original system - equals the scaled 
amount of work in an MjGjl queue in which the 
two customer classes are combined into one cus¬ 
tomer class. Based on standard heavy-traffic results 
for the M/G/l queue this implies that the distri¬ 
bution of the scaled total workload converges to an 
exponential distribution with mean pE[H^]/2E[i3], 
where E[H] = {pi+X 2 E[B'^])/{Xi+X 2 ) andEjH^] = 
(2Ai//rf + A 2 (Var[H'] + E[H']2))/(Ai + A 2 ). Since 
in HT almost all customers are located in Q 2 and 
^2 t ( 1 ~Pi)/E[-S 2 ]: the scaled number of customers 
in Q 2 is exponentially distributed with parameter 
r], 

PropertylSlis due to the the time-scale separation 
in HT, i.e., the dynamics of the stable queue evolve 
at a much faster time scale than the dynamics of 
the critically loaded queue. 


— + X 2 E[B' 2 ] = 1 - Slo* , w*>0,0<(5<l. 

Ml 

(14) 

Any arbitrary way in which we let Ai and A 2 ap¬ 
proach respectively A 4 and A 2 , for i5 j, 0 , will cause 
Q 2 to become critically loaded (because of assump¬ 
tion (TT^ l. All results obtained in this paper will 
still be valid, by choosing oj* = wEjH^]. 

In Section [3^ we argued that oj = l/EjH^] is a 
good choice, because it leads to dN 2 = (I — ■U 2 )A^ 2 - 
To see this, we first introduce the notation = 

(1 — pi)/E[B 2 ] and rewrite: 

limSN 2 = lim ( i(A"^*‘- A 2 ) | fVa- (15) 

HO A2tA^”‘ \0J ) 

I/E[i? 2 ] yields 

lim (-(Ar‘-A 2 ) ) fV 2 = lim ((1 - pi) - A 2 E[H'])) Af 2 

A2tA“"‘ \0J ) AatAf-* 

= lim (1 — U 2 ) N 2 , 

M2tl 

which has an appealing form. From (HU one can 
immediately see that oj = 1 results in taking the 
limit of (A 2 ’'**—A 2 )iV 2 , which is also quite commonly 
used (see for example [l^l. 

Finally, we consider increasing Ai and A 2 simulta¬ 
neously while keeping their ratio fixed. Denote by 
A = Ai -I- A 2 the total arrival rate, and let Aj = A^A 
(i = 1 , 2 ). As before, without loss of generality, we 
assume that Q 2 becomes unstable when increasing 
A. Taking oj = X 2 yields 


lim 

A2tA§"‘ 


(A 2 ="*‘ - A 2 ) N 2 = lim (A-^"* - K)N 2 


Substitution of w = 


4-2. Critical value and choice of oj 

In this note, we take the limit by increasing the 
arrival rate of Q 2 until it reaches its critical value 
while keeping the other parameters fixed. However, 
one could consider more general ways of varying the 
arrival rates in order to let Q 2 become critically 
loaded. To this end, we introduce A); and A2 such 
that Aj'/pi -I- A2E[i?2] = I. Additionally, we assume 
that 

fi2) 

ki l. + tjElB;]’ 

or equivalently: 



I 


(13) 


where A'^"* is the total arrival rate at which the 
second queue becomes critically loaded. 

4-3. Two critically loaded queues 

In the current note we have analysed the heavy- 
traffic behaviour under the condition Xi/ki < 
A 2 /fc 2 , i.e., only Q 2 becomes critically loaded. The 
limiting regime, in which both queues become sat¬ 
urated simultaneously (Ai/fci = X 2 /k 2 ), shows fun¬ 
damentally different system behaviour. In this case, 
it is not sufficient anymore to use a scaling that 
implies that exactly ^2 customers are served at Q 2 
during each cycle, i.e., the probability that there are 
served less than ^2 customers cannot be neglected, 
when analyzing the asymptotic behaviour of Qi. 
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5. Extensions 

The presented framework can be generalized in 
several areas without fundamentally complicating 
the analysis. 

5.1. General distributions 

Using the insights in the heavy-traffic system be¬ 
havior obtained from the perturbation analysis, a 
HT limit theorem for systems with general renewal 
arrival processes and general service and switch¬ 
over times can be conjectured. The interchange be¬ 
tween service and switch-over times proves its value 
again, since it makes it possible to link the workload 
of the system to the workload of a GIjGjl queue, 
for which the HT behavior is well-studied. Sub¬ 
sequently, the limiting distribution of the (scaled) 
number of customers in Q 2 can be derived. Let Ai, 
Bi, and Si denote the generic interarrival times, ser¬ 
vice times, and switch-over times for i = 1, 2, which 
can be generally distributed. Furthermore, let 

E[H'] = E[H 2 ] + ^, Var[B'] = Var[H 2 ]+^^^ 
k2 K2 

which is a generalisation to general service times of 
(g]), and let 

Ar‘ = (1 - Pi)/E[i?'] 

denote the limiting arrival intensity of customers at 
Q 2 . We scale the interarrival-time distribution A 2 
by changing its mean, while keeping the coefficient 
of variation fixed. We can formulate the following 
conjecture. 

Conjecture 1. The HT limit of the scaled queue 
length (1 — U 2 )N 2 as U 2 f 1 is exponentially dis¬ 
tributed with mean 

+ Ar'(VaHBa + c,4,l!:[S2l*)l), (16) 

where cv\. denotes the squared coefficient of varia¬ 
tion of Ai. 

5.2. N (> 2 ) queue polling system 

One could readily extend the analysis and results 
of the present note to an N{> 2) queue polling sys¬ 
tem with one critically loaded queue. In this case 
the stable queues have the same joint queue-length 
distribution as in an — 1 queue polling model 
with an extended switch-over time, while the scaled 
queue-length distribution of the critically loaded 
queue follows again an exponential distribution. 


6. Numerical results 

This section presents four numerical examples 
illustrating the validity and applicability of our 
proven limiting results and key concepts: 

1. (Joint) queue length distributions; 

2. Variance absorption result; 

3. Extension to general distributions; 

4. Extension to N{> 2) queue polling system. 

Example 1. In the first example we consider three 
polling models, referred to as PI, P2, and P3. In 
each model we have Ai = 0.06, pi = P 2 = 0.5, 
CTi = (72 = 0.5, ki = 2 and = 3. The only dif¬ 
ference between the three models are the values of 
A 2 , which are 0.25, 0.258, and 0.263, respectively. 
The corresponding values for U 2 are 0.953, 0.98, and 
0.997. Table [T] shows the simulated queue-length 
probabilities for Qi for these three models, and for 
the corresponding vacation model (model V). It 
can clearly be seen that the queue-length probabil¬ 
ities converge to those of the vacation system as 
stated in Theorem [T] 

Figure [ 1 ] shows how in the simulation the scaled 
queue-length (1 —^ 2)^2 in model PI closely follows 
the exponential distribution with parameter rj = 
1.289 as expected from Theorem!^ For models P2 
and P3 the differences between the simulation and 
the closed-form asymptotics are hardly noticeable, 
which is the reason why we have omitted them in 
the figure. 

Finally, we have calculated the correlation coef¬ 
ficient between the simulated number of customers 
in Qi and the scaled queue length of Q 2 . For PI, 
P2, and P3 the values are 0.048, 0.022, and 0.004, 
respectively. These values converge to 0, which 
follows from Property [3l implying that the amount 
of “memory” of the stable queue asymptotically 
vanishes compared to that of the critically loaded 


Model 

Po 

Pi 

P2 

P3 

Pi 

P5 

PI 

0.589 

0.275 

0.095 

0.029 

0.009 

0.003 

P2 

0.585 

0.277 

0.096 

0.030 

0.009 

0.003 

P3 

0.583 

0.278 

0.097 

0.030 

0.009 

0.003 

V 

0.582 

0.278 

0.097 

0.030 

0.009 

0.003 


Table 1: Simulated values for pj, := P[JVi = fc]. 
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Figure 1: The densities of (1 — U 2 )N 2 of model PI and an 
exponential distribution with parameter 77 . 


(1 - U 2 )N 2 

PI 

P2 

P3 

Mean 

0.742 

0.762 

0.778 

St. dev. 

0.757 

0.766 

0.781 

CV 

1.020 

1.006 

1.005 

(1 - p)N2 

PI 

52 

53 

Mean 

0.730 

0.757 

0.770 

St.dev. 

0.756 

0.768 

0.768 

CV 

1.035 

1.015 

0.996 


Table 2: Simulated means, standard deviations, and coeffi¬ 
cients of variation of(l — U 2 )N 2 for Example 2. 


queue. 

Example 2. As a by-product of our analysis we 
have observed that an interchange exists among the 
first two moments in service and switch-over times 
in systems with and without switch-over times such 
that the HT limits (for Q 2 ) are identical in both 
systems. That is, by noting that in HT every cycle 
precisely ^2 are served, the second terms in E[i? 2 ] 
and Var[i? 2 ] apportion the first two moments of the 
switch-over times to each of the customers fairly. It 
is important to note that this interchange is just 
one of the possibilities, i.e., a variety of trade-offs 
exists between service and switch-over times, all of 
which yield the same asymptotic distributions, as 
long as the terms remain unchanged. 

In this second example, we want to illustrate this 
generalization of the variance absorption result for 
branching systems Q with a numerical case. We 
start by considering again the three polling models 
PI, P2, and P3, but now we construct alternative 
models where (a portion of) the switch-over times 
is transferred to the service times of the customers 
in Q 2 . In more detail, we construct corresponding 
models with modified service times and switch-over 
times, satisfying 

1 l/tTl -I- 1/(72 

= — -h a - - -, 

^J -2 k 2 

1 , J/af + l/al 

2 T 7 7 

E[5i] = (1 - a)/(Ti, Var[5i] = (1 - a)/al 

E[52] = (1 - a)/(T2, Var[52] = (1 - a)/<jl 

choosing a = 1/2 for models PI, P2, P3, and o: = 1 
for models SI, S2, S3. Note that these latter three 
models are polling models without switch-over 




times and that M 2 = p in this case. The means, 
standard deviations, and squared coefficients of the 
scaled queue length {1 — U2)N2 are obtained using 
simulation, and given in Table [21 These values 
clearly confirm that the asymptotic distributions 
are identical, independent of the chosen inter¬ 
change, and that they converge to an exponential 
distribution with mean l/t] = 1/1.289 = 0.776. If 
the aforementioned interchange keeps the distribu¬ 
tion of the sum of ^2 service times in Q 2 and the 
switch-over times unaffected, then the HT limit of 
Qi remains unchanged as well. 

Example 3. In this example we use simulation 
to validate Conjecture [T] on the HT asymptotics 
for systems with general renewal arrival processes 
and general service and switch-over times. In 
many communication systems (some of these) 
distributions may be heavy-tailed, stressing the 
importance of knowing the limiting queue-lengths 
for general distributions. We take the same setting 
as in Example 1, but with different distribu¬ 
tions. Table [3] gives an overview of the various 
interarrival-time, service-time, and switch-over 
time distributions. Note that their means are 
equal to those in Example 1, where everything 
was assumed to be exponentially distributed. Also 
note that we have added a fourth system, P4, with 
A 2 = 0.2.635, for reasons that become clear later. 


Distribution 

Queue 1 

Queue 2 

Interarrival times 

Gamma(0.5,0.03) 

Pareto(2/(3A2), 3) 

Service times 

Uniform(0,4) 

Pareto(1.17, 2.4) 

Switch-over times 

Constant (2) 

Constant(2) 


Table 3: Settings for Example 3. 
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(1 — U 2 )N 2 

PI 

P2 

P3 

P4 

Mean 

7.5 

18.4 

116.3 

232.5 

St.dev. 

30.5 

38.6 

192.1 

281.8 

CV 

4.0 

2.1 

1.7 

1.2 


Table 4: Simulated means, standard deviations, and coeffi¬ 
cients of variation for Example 3. 


We have simulated the queue-length distributions 
of each of these four systems and compare them to 
the limiting values, based on the exponential distri¬ 
bution that arises in the HT limit. The parameter 
of this exponential distribution, which can be com¬ 
puted using (ITC)) . is 77 ~ 2.527. Table S] gives the 
mean, standard deviations, and coefficients of vari¬ 
ation of the scaled queue lengths for each of the 
models. From this table we conclude that the first 
two moments indeed approach those of an exponen¬ 
tial distribution, albeit in a slower pace than Ex¬ 
ample 1, which had only exponential distributions. 
This slower speed of convergence is the reason for 
including a fourth model, P4, with an even higher 
utilization. Figures [ll[a)-(b) shows the simulated 
scaled queue length probabilities for PI, ..., P4, 
and the density of the limiting distributions, again 
confirming the convergence to the exponential dis¬ 
tribution. 





Figure 2: Simulated and theoretical scaled queue lengths for 
Example 3. 


Example 4. In this final example we illustrate 
how queue lengths behave in a polling system with 
N{> 2) queues, each with fc-limited service, when 
all arrival rates are increased simultaneously, as dis¬ 
cussed in Section O We consider a system with 
N = A queues and exponentially distributed ser¬ 


vice times, interarrival times, and switch-over times 
with /ii = 1/2, ffi = 1/3, for i = 1,2,3,4. The ser¬ 
vice limits are ki = 7, k 2 = ks = 6 , and k^ = 5. Let 
A = Ai • -I- A 4 denote the total arrival rate. We 
will increase A until, one by one, all of the queues 
become unstable. While increasing A, we keep the 
ratio between the four arrival rates fixed. Let Xi 
denote the proportion of the total arrivals that is 
routed to Qi, i.e. 

A. = A,A, (17) 


with Ai = i/10 in this example. The impact of 
increasing the total arrival rate on the system be¬ 
havior is graphically illustrated in Figure [31 which 
depicts simulated mean queue lengths E[W(A)] as 
a function of the total arrival rate A. 



Figure 3: The mean queue lengths in Example 4 as functions 
of the total arrival rate A. The four vertical dashed lines 
correspond to respectively A|"‘, Ag"^, AJ"* (from left 
to right). 


The following interesting observations can be made: 

• Obviously, for A = 0 all mean queue lengths 
are 0 . 

• As A is increased, all mean queue lengths in¬ 
crease. Queue 1 has the lowest utilization, 
causing the mean queue length to increase 
slowest. Queue 4 has the highest utilization 
and, as a consequence, is the first queue to 
become unstable, which happens at A 4 "* = 
25/74 Ri 0.34 (the leftmost dashed vertical 
line). This value can be computed via (HI and 

m- 

• Note that, although Q 4 is unstable at A = 
A 4 "*, the other queues are stable. We can 
increase A further, causing the mean queue 



































lengths of the other queues to increase, while 
(54 remains unstable. 

• When A approaches the value of A 3 ’'®* = 
10/23 « 0.43, Qz becomes unstable as well. 
Now Qi and Q 2 are the only stable queues. 

• As we continue to increase the total arrival 
rate, Q 2 and, eventually, also Qi become un¬ 
stable. This happens at respectively A 2 ’'®* = 
15/26 Ri 0.58 and AJ®-®* = 7/6 ss 1.17. 

• The mean queue lengths (as functions of A) are 
not differentiable at A = {AJ’'®*,..., A|’'®*}. 

This particular behavior can be explained using the 
results from this paper. As discussed in Section 
15.21 from the viewpoint of customers in queues 1, 
2, and 3, the four-queue polling system behaves like 
a three-queue polling system when (54 is unstable, 
i.e. for Al®"®* < A < Ag’'®*. In this so-called “cor¬ 
responding system” the switch-over time between 
Qs and Qi is the sum of the switch-over times 
S'3 -I- 5'4 in the original model, plus service times 
that are exp(/i 4 ) distributed. In this three-queue 
polling model, the next queue to become unstable 
is ( 53 , at Ag®"®* = 10/23 ^ 0.43, making the sys¬ 
tem behave like a two-queue polling model with 
an even longer switch-over time between Q 2 and 
Qi. This continues, with Q 2 becoming unstable at 
A = 15/26 Ki 0.58 and, finally, Qi being the last 
queue to become unstable when the total arrival 
rate reaches the value of AJ®"®* = 7/6 « 1.17. 

The simulation results can be used to verify our 
claims that, near the point where one of the queues 
becomes unstable, (a) the number of customers in 
the stable queues are distributed as in a system 
with that particular queue removed and replaced 
by an extra long switch-over time, and (b) the 
scaled queue length of the unstable queue is ex¬ 
ponentially distributed. Therefore, we consider the 
original polling system near each of the critical val¬ 
ues, namely at A = 0.995A/®'®*, for all i, and com¬ 
pare the simulation results to the simulation results 
of the “corresponding” polling systems. Table [5Ka) 
shows the simulated queue length probabilities of 
the stable queues. For A A 4 ®"®* queue 4 becomes 
unstable and the system starts behaving as a three- 
queue polling system with an extra long switch-over 
time S 3 . The simulated queue length probabilities 
of this “corresponding system” are listed in Table 
[3(b). It is readily seen that the results are com¬ 
pletely equivalent, the extremely minor differences 


being attributed to simulation inaccuracies. This 
has been repeated near A r: Ag®'®* and A r AI®'®*, 
and the same conclusion can be drawn. 





Ar 

A|®'®‘ 



Queue 

Po 

Pi 

P2 

P3 

P4 

P5 

<5i 

0.534 

0.292 

0.118 

0.040 

0.012 

0.003 

Q2 

0.325 

0.281 

0.186 

0.106 

0.055 

0.026 

Qs 

0.186 

0.205 

0.177 

0.136 

0.098 

0.067 





Ar 

A crit 



Queue 

Po 

Pi 

P2 

P3 

P4 

P5 

<5i 

0.410 

0.301 

0.166 

0.076 

0.031 

0.011 

Q2 

0.177 

0.200 

0.179 

0.142 

0.103 

0.071 





Ar 

Ay®‘ 



Queue 

Po 

Pi 

P2 

P3 

P4 

P5 

<5i 

0.296 

0.273 

0.195 

0.119 

0.064 

0.031 


(a) Original polling system 



Table 5: Marginal queue lengths of the stable queues in the 
original polling and the corresponding systems. The simu¬ 
lated probabilities ¥(Ni = m), denoted by Pm, are given 
for m = 0,1, 2,. .., 5. 

Verifying that the scaled queue lengths of the un¬ 
stable queues are exponentially distributed, which 
is our second claim, requires determining the pa¬ 
rameters of these distributions as a first step. These 
parameters can be computed similar to (ITBl) . but 
with i ?2 modified to include the service times and 
switch-over times of all overloaded queues, and Bi 
the service time of an arbitrary customer in any of 
the stable queues. The values of rn, i = 1,2, 3 , 4 , 
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equal 1.364, 2.875, 3.739 and 4.529 taking oj = Xi. 
The latter is the appropriate value when consider¬ 
ing the scaling — h)Ni as suggested in Section 
14.21 Figures l4l[a)-(d) depict the simulated proba¬ 
bility density functions of the scaled queue lengths 
(Acr*t_ A) 7 Vi, for i = 1,... ,4 at A = 0.995Af’'*‘, i.e. 
very close to instability. It is clear that the simu¬ 
lated values are extremely close to the theoretical 
values. 



Figure 4: Simulated and theoretical scaled queue lengths for 
Example 4. 

To obtain more insight in the influence of over- 
loaded queues, we have chosen to do focus on 
one specific part of Figure |31 namely the part be¬ 
tween A = Ag''** and A = A 2 ’'**, i.e. the interval 
[0.435,0.577]. In this interval, queues 3 and 4 are 
unstable, causing the system to behave like a two- 
queue polling system from the perspective of cus¬ 
tomers in queues I and 2. This corresponding sys¬ 
tem is similar to the (first two queues of the) origi¬ 
nal system, but has a longer switch-over time from 
Q 2 to Qi, which is distributed as the convolution of 
the original switch-over times S 2 , 83 , 84 , and fcg ser¬ 
vices of type-3 customers, and ^4 services of type-4 
customers. Figure [S] depicts the mean queue lengths 
of queues I and 2 in the original model (taken from 
Figure El) and, additionally, in the corresponding 
two-queue polling model, while varying the total 
arrival rate A between 0 < A < 0.577. It can 
clearly be seen in Figures EKa) and EKb) that the 
mean queue lengths of the two systems are differ¬ 
ent for A < 0.435, but they are the same in the 
aforementioned region 0.435 < A < 0.577. This 
confirms that the two systems behave identically in 
this region. 


E[A'l(A)] Ag"‘ A§"‘ 




Figure 5: Simulated mean queue lengths for queues 1 and 2 
in the original model (solid lines) and in the “corresponding” 
polling model (dashed lines). 

7. Conclusion 

The present note has studied a novel applica¬ 
tion of the singular-perturbation technique leading 
to a rigorous proof of the HT limits in A:-limited 
polling models with switch-over times. The scaled 
queue-length of the critically loaded queue follows 
an exponentially distribution of which the parame¬ 
ter is known. The number of customers in the sta¬ 
ble queue has the same distribution as the number 
of customers in a vacation system. Vacation sys¬ 
tems with fc-limited service can be analyzed very 
efficiently using modern numerical techniques. In 
particular, when the service times and switch-over 
times are constant, or in case they have a phase- 
type distribution or a Gamma distribution, the no¬ 
torious root-finding problem becomes trivial (see 
for more details). As such, the asymptotics 
form an excellent basis for approximating fc-limited 
polling models - and, thus, WPAN and mobile (ad- 
hoc) networks - with general load. 
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